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Dynamical 2+1 flavour lattice QCD is used to calculate the splittings between the masses of 
mesons and baryons containing a single static heavy quark and domain-wall light and strange 
quarks. Our calculations are based on the dynamical domain-wall gauge field configurations 
generated by the RBC and UKQCD collaborations at a spatial volume of (2.7 fm) 3 and a 
range of quark masses with a lightest value corresponding to a (partially-quenched) pion mass 
of 275 MeV. When extrapolated to the physical values of the light quark masses, the results 
of our calculations are generally in good agreement with experimental determinations in the 
bottom sector. However, the static limit splittings between the f2^~ baryon and other bottom 
hadrons tend to slightly underestimate those obtained using the recent D0 measurement of 
the 



I. INTRODUCTION 



Hadrons containing a single bottom quark have received much attention recently. Over the last 
decade, the bottom meson sector has been investigated in great detail at the b factories (BELLE and 
BaBar) and using the TeVatron at Fermilab. With the recent D0 measurement of the baryon 
p], we are rapidly approaching a complete picture of the ground-state bottom baryon sector as 
well. The start up of the Large Hadron Collider (LHC) will dramatically increase our knowledge 
of b physics; both the dedicated b physics experiment, LHCb, and the general purpose ATLAS 
and CMS experiments expect to observe unprecedented numbers of bb pairs, many of which will 
hadronise to bottom mesons and baryons. 

In the bottom meson sector, many important observations have been made in the last decade, 
dramatically refining our knowledge of flavour physics. The CKM mechanism of the Standard 
Model currently provides a good description of current flavour-changing measurements of B and 
B s mesons, bringing us into the precision era of flavour physics [2]. The various LHC experiments 
aim to investigate the decays of bottom baryons with enough precision to further test the Standard 
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Model. In particular, the fact that bottom baryon polarisation can be measured may uncover new 
right-handed couplings. To search for physics beyond the Standard Model in these systems, it 
is necessary to have accurate predictions from the Standard Model. Typically, this requires non- 
perturbative evaluations of matrix elements from lattice QCD and thus an understanding of bottom 
hadrons in lattice QCD. 

The simplest properties of the bottom hadrons are their masses. Before more complicated bot- 
tom observables can be predicted with any rigour, it is necessary to have reliable lattice calculations 
of the masses and mass splittings. There are many lattice studies of the mass of the B mesons in 
the literature but fewer of the bottom baryons [3] H] El E] and, only recently, the first unquenched 
bottom baryon calculations have appeared [10]. In this paper, we use light up and down 
quark masses to study the spectrum of bottom baryons and mesons. For the bottom quark, we 
work in the static (m& — > oo) limit. Since the bottom quark mass, m;, 3> Aqcdj many properties 
of the physical bottom hadrons are expected to be close to those of hadrons containing a single, 
infinitely massive (static) quark, with corrections suppressed by powers of Aqcd/w-6- Such correc- 
tions can be addressed systematically in Heavy Quark Effective Theory (see, e.g. Ref. [H] for a 
review). Our calculations use domain- wall fermions for the light and strange quarks and make use 
of the dynamical gauge configurations generated by the RBC and UKQCD collaborations using 
domain-wall quarks |12| [13]. We currently work at a single lattice spacing, a = 0.114 fm |13] and 
focus on a volume of spatial side, L = 2.74 fm, for a range of quark masses corresponding to pion 
masses between 275 and 750 MeV. 

This paper is organised in the following manner. In Section In] we present the details of the 



lattice calculations we have carried out. We discuss the mass splittings we observe in Section III 



before concluding in Section IV 



II. DETAILS OF LATTICE CALCULATIONS 



The lattice QCD calculations presented here are based on the ensembles of 2+1 flavour lattices 
generated by the RBC and UKQCD collaborations |12[ [13]. These ensembles use the Iwasaki 
gauge action |14[ [T5] with (5 = 2.13 and a domain wall quark action using the length of the fifth 
dimension, L s = 16. The parameters of the ensembles used this work are shown in Table [I] and we 
refer the reader to Ref. |12[ [T3] for further details. 

Using these ensembles, we have computed domain wall quark propagators |16| [TT] ITS] ITU] EO] 
for various different (partially quenched) quark masses as shown in Table Ull Our inversions use 
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Ensemble 




Volume 


(sea) 

TO ii g ht 


(sea) 
^strange 


m res [T3] 


A cfgs 


A 


2.13 


24 3 x 64 


0.005 


0.04 


0.00315 


- 140 


B 


2.13 


24 3 x 64 


0.01 


0.04 


0.00315 


~ 180 


C 


2.13 


24 3 x 64 


0.02 


0.04 


0.00320 


- 125 



TABLE I: Parameters of the ensembles of RBC/UKQCD gauge configurations used in this calculation. For 
full details, see the original works, Refs. [12| [13]. 



L s = 16 and a domain wall height of M5 = 1.8 (note that these are the same parameters used in 
generating the ensembles and the bold entries in Table [TT] correspond to QCD computations). To 
obtain clean signals, for the hadron energies and splittings, we use APE [21, 22 smeared sources 
at multiple locations on each gauge configuration, performing a separate inversion for each source. 
The approximate masses of the pseudoscalar mesons computed with various combinations of these 



propagators are shown in Table III 



Ensemble 


(val) 
TO light 


^cfg 


N SIC 


A 


0.002 


- 140 


5 


A 


0.005 


- 140 


6 


A 


0.01 


- 140 


1 


A 


0.02 


- 140 


1 


A 


0.03 


- 140 


1 


A 


0.04 


-140 


6 


B 


0.005 


- 180 


1 


B 


0.01 


-180 


5 


B 


0.02 


- 180 


1 


B 


0.03 


- 180 


1 


B 


0.04 


-180 


5 


C 


0.02 


-125 


1 


C 


0.04 


-125 


1 



TABLE II: Quark propagators used in calculations of the mass splittings. Bold entries denote rn^^J — 
m Hght 1 • -^sre indicates the number of sources that were used on each gauge configuration. These sources 
were spread around the lattice volume at as large relative separations as possible. 

The bottom quark is implemented in the static limit. Its propagator, Sq is represented as a 
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Ensemble 


(val) 
m light 


(val) 
^"strange 


mi val) [GeV] 


m£ al) [GeV] 


A 


0.002 


0.04 


0.275 


0.560 


A 


0.005 


0.04 


0.331 


0.576 


A 


0.01 


0.04 


0.415 


0.602 


A 


0.02 


0.04 


0.546 


0.654 


A 


0.03 


0.04 


0.653 


0.701 


A 


0.04 


0.04 


0.747 


— 


B 


0.005 


0.04 


0.335 


0.581 


B 


0.01 


0.04 


0.419 


0.607 


B 


0.02 


0.04 


0.550 


0.657 


B 


0.03 


0.04 


0.657 


0.706 


B 


0.04 


0.04 


0.751 




C 


0.02 


0.04 


0.549 


0.654 



TABLE III: Approximate pion and kaon masses for each set of measurements. 



product of gauge links in the temporal direction 

S Q {x,t;t ) 1 



t'=t 



(1) 



where are SU(3) gauge links. At non-zero lattice spacing, there are different discretisations of 
the heavy quark action that have the same continuum limit. It is well-known (see e.g. Ref. [23 ) 
that signals for static hadron quantities are improved if the gauge links appearing in Eq. ([!]) are 
smeared in some manner over a small local volume. After extensive testing to optimise the heavy 
hadron signals, we perform our calculations with hypercubically-smeared (HYP) gauge links |24j 
in the heavy quark propagator, Eq. ([I]). We study a number of choices of heavy quark smearing 



parameters as shown in Table IV labelled Si for i = 0, . . . , 5. 



Set 


So 


Si 


S 2 


s 3 


S4. 


s 5 


Ahyp 


10 


5 


1 





3 


3 


ai 


0.75 


0.75 


0.75 


0.75 


0.75 


0.6 




0.75 


0.75 


0.75 


0.75 


0.4 


0.6 


a 3 


0.75 


0.75 


0.75 


0.75 


0.65 


0.6 



TABLE IV: Parameters used in HYP-smearing the heavy quark action. 
To extract the lattice energies, Eh, of the various hadrons, h, we compute the meson and baryon 
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two-point correlation functions 

C/(Mo) = £tr[s Q (x,t;io)S}(x 5 f;x»<o)] , (2) 

X 

C f, 9 & ^) = *O x ' <5 *o)<*' V'"*' (sf (x, t; x, i )r) (rsf (x, t; X, *„)) , (3) 



for the various light flavour combinations, f(g), and then combine them appropriately. The trace 
in the meson correlator is over color and spinor indices, and where explicit, the upper(additional 
lower) indices on the propagators correspond to colour (spin). For the baryons with the spin of 
the light degrees of freedom, S£, being zero we choose T = C75 where C is the charge conjugation 
matrix, while for the baryons with s# = 1, we measure the three polarisations, corresponding to 
r = C7123, and average them in extracting the energies as they are degenerate in the infinite 
statistics limit. We combine the various correlators in Eq. ([3]) as appropriate for the specific SU(3) 
representation. For example, the belongs to the flavour anti-triplet and is a si = state, 
consequently 

C S o = CZ% - C]% (4) 

and similarly, 

Cao = CT» - Cj^=0^ + (%, CUo = C& + C2» B , C =C]y (5) 

' ' 0,1 > J b,i 0,1 

As we work in the isospin limit, the other singly- heavy hadrons are degenerate with those above. 
In the static limit of the heavy quark, the J p = ^ + states E' b , Q&) and | + states (E^, H^*, fi^) 
are degenerate and our results correspond to this. Future calculations for 7^ 00 will enable the 
spin splittings to be investigated, see Ref. [8.1 for recent work. 

The transfer matrix formalism immediately shows that in the limit of large temporal separation, 
these correlators are dominated by exponentially decaying signal of the ground-state hadron energy. 
That is, 

C h it, t ) ^ A h exp [-E h (t - t )} , (6) 

where is the lattice energy. In the case of the S£ = 1 baryons of strangeness < 2, this is 
complicated by the fact that, at the physical quark masses and in infinite volume, they decay 
strongly to lighter s/> = 0,1 baryons and a pion. For example, — ► A(,7r + is the dominant 
decay mode [25 . However, at the unphysical masses and finite volumes used (the pion couples 
derivatively, requiring non-zero lattice momentum) here, these decays are forbidden. Additionally, 
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the finite temporal extent (non-zero temperature) of the lattice geometry used in our calculations 
results in the appearance of exponentially suppressed states that decay with "lighter masses" that 
may pollute the signals of the zero temperature ground-state [26\ (this is asymmetric analog of the 
constant term expected in a two meson correlation function from one meson propagating forward 
in time and one propagating backward in time |27l EE])- We will return to this issue below. 



III. BOTTOM HADRON MASS SPLITTINGS 



In the static limit, the masses of bottom hadrons are infinite, however, the residual mass, A^, 
of a heavy hadron h is well defined. In the lattice regularisation used here, the lattice energies, 
Eh diverge as the continuum limit is taken and a perturbative subtraction is required to convert 
them to the corresponding residual mass, for example, in the modified minimal subtraction (MS) 
scheme. At one loop, this can be performed straightforwardly for typical lattice actions, and with 
more effort for the HYP-smeared gauge links used herein |29[ [30] . However, there are significant 
issues with higher order corrections and renormalons in the subtraction procedure |31[ [32] . In 
contrast, the differences between the residual masses of two hadrons containing the same heavy 
quark (that is, the same lattice discretisation of a static quark) have well-defined continuum limits 
in which they are observable (the perturbative corrections and renormalisation scale dependence 
cancel in the difference). In our work, we concentrate on these differences and thereby avoid the 
issues mentioned above. 



A. Heavy quark action 

We first study the effects of the choice of the heavy quark action via the application of the 



various HYP-smearings shown in Table IV For each energy or energy difference, the results for 



all the choices of smearing are found to be consistent within uncertainties. Figure [T] shows the 
noise-to-signal ratios for the correlators corresponding to the — B u , — B s and — A?, mass 
splittings for the various smearing parameter sets as a function of Euclidean time using the QCD 
propagators of ensemble B in Table [IJ As can be seen, the parameter set Sq clearly has the strongest 
signals. The other sets of propagators display similar behaviour and we focus on set Sq for the 
remainder of our study. 
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FIG. 1: The ratio of noise to signal, Rns, for exemplary correlator differences for the different choices of 



the heavy quark action (Table IV I. The m( val ) = m( sca ) quark propagators on the B ensemble were used in 



this comparison. The other sets display similar behaviour. 

B. Analysis Methods 

To extract the mass differences from the lattice correlators defined in the previous section, 
we perform multiple independent analyses using different methods. We fit to either single- or 
two-state effective masses and /or mass differences (see Ref. [33]), or perform multi-exponential 
fits to (differences of) correlation functions using Bayesian priors (as in Refs |344 [3l)]). We use 
the bootstrap analysis method to estimate our uncertainties. We average over all propagators 
computed from the various sources on each configuration and then bin the measurements over 
units of ~100 molecular dynamics time units as residual correlations are observed to persist over 
this range [13] . The number of bootstrap ensembles used in our analysis is typically four times the 
number of binned measurements, bootstrap ~ 200. 

Correlation functions are calculated with both an APE smeared sink (using the same parameters 
as at the source) or a point sink. Results from both sink types are consistent. In general, the fitting 
to the energy differences leads to a cleaner signal than fitting each energy individually, but not in 
all cases and both approaches were tried. 

In the correlated fits to single- and two-state effective masses, a systematic uncertainty is as- 
signed from varying the fitting ranges, imin to ^max? or. from differences between fits to sliding 
windows of time-slices within the overall fit range. Typically, single effective masses are fit from 
^min ~ 9 and two-state effective masses are fit from t m j n ~ 4 (here, the ground state plateaus 
earlier as the excited state contamination is decreased significantly). The upper limit of the fits 
is imax ~ 20 for mesons and anti-triplet baryons, and i max ~ 15 for sextet baryons where states 
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arising from the finite temporal extent discussed above pollute the signal. 
In the Bayesian analysis, the correlation functions are fit to 

C(t) = A 1 e- E ^(l + B 1 e~ AE ' t + B 2 e~^ +AE ^ t )+A 1 C 1 e~ E ' t . (7) 

over the range 2 < t < 20. This form was sufficient to provide an acceptable fit to all correlators and 
additional exponential terms were not constrained by the data. The Bayesian prior distribution 
functions are Gaussian. The mean values for the A\ and E\ priors were estimated from single 
exponential fits at large times. The widths for those parameters were taken to be about 50% and 
10%, respectively, are significantly larger than the uncertainties from the initial single exponential 
fit. The mean values for the B{ priors were in [—1, 1] with widths up to a factor 10, and the AEi 
priors had means in [0.3,0.6] with widths between 60% and 80%. The last term in allows 
for contributions to C{t) due to multi-particle states created by the interpolating operators in 
Eq ([2]) and Q. Usually these states give vanishingly small contribution, as their total energy is 
greater than the hadron mass of interest. However, with (anti-) periodic boundary conditions in t, 
these states can contaminate the large t behavior due to a backward propagating light meson [26 . 
Therefore we include the last term in Eq. ([7| with prior mean for C\ ~ exp(— 0.1T), where T is 
the temporal extent of the lattice, and for E' = 0.5, both with 100% widths. 

In combining the various analyses, we make a direct comparison of the various extractions and 
if they are in agreement, we take their mean. If there are discrepancies (very few are found), we 
assign an additional systematic uncertainty large enough to make the result consistent with the 
individual analyses. 

The final results of our analyses of the lattice calculations are shown in Fig. [2] for the various 
quark masses on the three ensembles (the statistical, and all systematic uncertainties have been 
combined in quadrature). The empty points correspond to calculations that are partially quenched, 
with the sea quark masses differing from the valence quark masses. Ensembles A, B and C are 
coloured blue, red and green respectively. 

C. Quark mass extrapolations 

The light quark mass dependence of the static hadron energies and mass differences calculated in 
this work is, at least in principle, described by the low energy effective theory, heavy hadron chiral 
perturbation theory (HHxPT) [361 EH ESI ESI SO] or its quenched and partially-quenched versions 
|41l W2\ H3] . In addition to polynomial dependence on the pion and kaon masses, non- analytic terms 
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such as log[m^//u 2 ] appear (/i is the renormalisation scale). We have attempted to perform 
extrapolations using the forms predicted by these theories. Unfortunately, our current data are 
insufficient to constrain the coefficients of the non-analytic terms in the HHxPT expressions (three 
axial couplings contribute in the coefficients of the non-analytic terms in the various differences). 
Future, direct calculations of the axial couplings will enable a more controlled chiral extrapolation. 

With this in mind, here we use simple polynomial fits to perform the light quark mass extrapola- 
tions, allowing linear, quadratic and cubic dependence on m q ~ m 2 . Since we keep m s fixed, we do 
not include dependence on this parameter (we note that it is not quite tuned to the physical value 
|13|). These fits use a maximum likelihood estimator and are performed in an uncorrelated manner 
for simplicity. Accounting for the correlations between the different partially-quenched data points 
computed on the same underlying ensemble would slightly reduce our uncertainty. However, given 
the extrapolation forms used at present are ad hoc, we do not pursue this further. The resulting 
X 2 s of the fits are all acceptable. 

We have also performed coupled fits to all the differences, allowing a fixed polynomial depen- 
dence on ttItt for each of the hadron energies. This coupled fitting procedure provides a successful 
fit and agrees with the above (uncoupled) method within uncertainties but does not improve the 
extrapolations. 

In Fig. [2] we show the results of the linear (blue), quadratic (red) and cubic (green) fits to the 
data contained in the ordinate extent of corresponding shaded region. The extracted splittings at 
the physical light quark masses are presented in Table [V] for each extrapolation and compared to 
experimental determinations where available. Differences between the three forms of extrapolations 
are relatively small, and, as an example, the linearly extrapolated mass differences are shown in 
Fig [3] along with the experimental determinations. 

IV. CONCLUSIONS AND OUTLOOK 

In this paper, we have presented a study of the mass splittings between mesons and baryons con- 
taining one static quark using lattice calculations based on the domain- wall fermion RBC/UKQCD 
ensembles. We have used domain wall quarks with a range of masses, the lightest set corresponding 
to a pion mass of ~275 MeV and at a single lattice spacing, a = 0.114 fm. In general, the splittings 
we extract agree well with experiment. Those involving the O^" tend to be smaller than found 
experimentally as also found recently in the lattice calculations of Refs. [HI [§], but the discrep- 
ancy is barely significant at our current level of precision once the quark mass extrapolations are 
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hi 


h 2 


M hl - M h2 (lin.) 


M hl - M h2 (quad.) 


M hl - M h2 (cub.) 


M hl - M h2 (expt.) 


B s 


B u 


0.097(28) 


0.092(35) 


0.084(32) 


0.0867(11) 


A b 


B„ 


0.362(88) 


0.37(11) 


0.36(12) 


0.3407(20) 


A 6 


B a 


0.271(88) 


0.28(12) 


0.29(10) 


0.2540(22) 




B u 


0.579(69) 


0.595(94) 


0.588(98) 


0.5136(30) 




B s 


0.482(55) 


0.500(75) 


0.502(82) 


0.4269(32) 




A,, 


0.214(67) 


0.211(86) 


0.212(76) 


0.1729(36) 


y+ 


B u 


0.600(68) 


0.564(89) 


0.556(93) 


0.5322(30) 


T+ 
^b 


B s 


0.508(64) 


0.484(88) 


0.491(96) 


0.4455(32) 


^b 


A,, 


0.222(90) 


0.19(13) 


0.21(14) 


0.1915(36) 


y+ 
^b 


-0 
^6 


0.012(80) 


-0.02(11) 


-0.002(93) 


0.0186(42) 


"6 


B u 


0.718(71) 


0.698(98) 


0.677(90) 


— 


77 '0 
"6 


B s 


0.633(64) 


0.622(91) 


0.616(79) 


— 


77>0 
"6 


A b 


0.332(94) 


0.30(13) 


0.27(13) 


— 


77'0 
"6 


77 o 


0.131(78) 


0.10(11) 


0.094(95) 


— 


"'0 
"6 




0.102(50) 


0.094(64) 


0.079(92) 


— 


fi- 


B„ 


0.817(42) 


0.794(61) 


0.776(58) 


0.886(15) 




B s 


0.718(37) 


0.705(48) 


0.696(51) 


0.799(15) 


n 6 - 


Ab 


0.419(78) 


0.385(98) 


0.37(11) 


0.545(15) 


n 6 - 


77 


0.222(65) 


0.182(97) 


0.175(83) 


0.372(15) 


fi- 




0.196(61) 


0.181(82) 


0.173(55) 


0.354(15) 


n 6 - 


77/0 
"6 


0.086(43) 


0.084(47) 


0.076(55) 





TABLE V: Extracted static limit mass differences [GeV] obtained through linear, quadratic or cubic ex- 
trapolations in TO^ to its physical value. The uncertainties in the extrapolated values correspond to one 
standard deviation allowed region of the space of fit parameters. Experimental mass splittings [TJ [55] are 
shown for comparison (for the mesons and the sextet baryons, we report the lighter of the heavy quark spin 
multiplcts). The SJ, has not been observed. 

accounted for. 

Our chiral extrapolations are somewhat ad hoc at present as the axial couplings between b- 
hadrons and light pseudoscalar mesons that appear in HHxPT are relatively unknown, making 
extrapolations using the results of chiral perturbation theory problematic. Ongoing calculations of 
these charges will allow us to control better the mass difference extrapolations. Additionally, the 
effects of the slightly unphysical strange quark mass need to be accounted for. Our calculations are 
performed in the static limit of the b quark and so our extracted values for the splittings are subject 
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FIG. 2: Polynomial extrapolations of the mass splittings calculated on the various ensembles are shown as a 
function of the pion mass. The blue, red and green points correspond to ensembles A, B and C respectively. 
Solid (open) symbols denote QCD (partially-quenched QCD) calculations. The uncertainties on the lattice 
data combine the statistical uncertainties and the various systematic uncertainties in quadrature. The blue, 
green and red shaded regions show the one standard deviation allowed region of linear, quadratic and cubic 
fits to the data within the ordinate extent of the corresponding region. 
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FIG. 3: SU(3) mass splittings between bottom hadrons linearly extrapolated to the physical light quark 
masses. The green shaded regions correspond to the experimental determinations (combinations with both 
the B* s and the are shown in purple). 

to uncertainties 0(AQCD/ m f>) when compared to the experimental bottom hadron spectrum. In 
the future, we intend to revisit these calculations using either non-relativistic QCD or the Fermilab 
heavy quark formalism. We also intend to extend our calculations to other lattice spacings and 
volumes as they become available. 
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